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Abstract. We consider dark matter in a minimal extension of the Standard Model (SM) 
which breaks electroweak symmetry dynamically and leads to a complete unification of the 
SM and technicolor coupling constants. The unification scale is determined to be M\j ~ 2.2 x 
10^^ GeV and the unihed coupling au « 0.0304. Moreover, unihcation strongly suggest that 
the technicolor sector of the model must become strong at the scale of O(TeV). The model also 
contains a tightly constrained sector of mixing neutral fields stabilized by a discrete symmetry. 
We find the lightest of these states can be DM with a mass in the range ut-dm ~ 30 — 800 GeV. 
We find a large set of parameters that satisfy all available constraints from colliders and from 
dark matter search experiments. However, most of the available parameter space is within 
the reach of the next generation of DM search experiments. The model is also sensitive to a 
modest improvement in the measurement of the precision electroweak parameters. 

Keywords: Dark matter, Naturality, Unification, Technicolor 


Preprint: HIP-2015-12/TH 


Contents 


1 Introduction 1 

2 General model of dark matter 3 

2.1 Gauge interactions 3 

2.2 Scalar interactions and masses 4 

2.3 Rotation to the mass eigenbasis 5 

2.4 Renormalizable and anomaly free implementation 6 

3 Unification 7 

3.1 Unification of the SM coupling constants 7 

3.2 Unification of all couplings including 04 = otc 9 

4 Experimental constraints on low energy theory 11 

4.1 Results of a generic MCMC scan 14 

4.2 Benchmark models 16 

5 Additional bounds: Higgs to 77 and a light composite scalar 17 

5.1 Higgs to 77 18 

5.2 Dynamical Higgs boson mass scale 19 

6 Conclusions 21 


1 Introduction 

In contrast with its celebrated finding of the Higgs boson [1, 2], LHC has so far failed to 
discover any obvious sign for new physics. Yet compelling evidence for new physics abound. 
For example the existence of a large dark matter component in the energy density of the 
universe is currently lacking a proper elementary particle physics context. Standard model 
is also plagued by the hierarchy problem, i.e. it lacks a natural explanation for the lightness 
of the Higgs particle. Finally, SM does not give rise to coupling constant unification, which 
would be highly desirable on theoretical grounds. Yet another issue is that Standard Model 
cannot explain why the universe contains only matter. Of these issues the DM problem is 
the most tangible one phenomenologically, and DM studies have indeed recently been the 
leading motivation for the construction and analyses of beyond SM scenarios. 

Dark matter can of course be considered separately from the issues of hierarchy problem 
and unification. In particular scalar dark matter models have been very popular recently [3- 
10]. These models are attractive because of their simplicity and because of possible connection 
to a new, larger dark sector interacting with the SM through the Higgs portal. Alternatively 
one can require that the DM is a part of a larger, independently motivated particle physics 
model. In this class the leading paradigm until recently was supersymmetry (SUSY) and 
in particular the minimal supersymmetric standard model. SUSY is compelling for many 
reasons, which include theoretical connections to (super)string theory and supergravity. From 
the phenomenological point of view SUSY is attractive because it could provide a solution 
to all main issues mentioned above. However, no sign of SUSY has been found. Instead 
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the minimal supersymmetric standard model is getting ever more tightly constrained by the 
data. 

Another popular model building paradigm involves new strong dynamics (technicolor) 
sourcing the electroweak symmetry breaking. The currently favored technicolor (TC) models 
are based on the idea of quasiconformality [11, 12], concrete realizations of which are the 
minimal and next to minimal walking technicolor models [13, 14]. Technicolor was originally 
put forward as a solution of the hierarchy problem. However, it was recently shown that TC 
models can easily also provide unihcation of SM gauge couplings [15-17] and they have also 
shown to contain several possible DM candidates [15, 16, 18-29]. 

In this paper we consider a model for particle DM, which is deeply motivated by the 
TC paradigm and by the requirement of the gauge coupling unification. Yet, the DM-sector 
of the model may be considered also independent of its TC context. Therefore we first 
discuss the generic low-energy setup featuring new leptons transforming as a doublet and a 
triplet under the weak gauge group, and their mixing patterns and couplings with the gauge 
fields and with a Higgs-like scalar sector. The setup naturally leads to a thermal relic DM 
candidate, completely external to the strongly interaction sector and with fully deterministic 
couplings with ordinary matter. We also discuss a concrete model for dynamical electroweak 
symmetry breaking, which determines all mass terms and effective Higgs couplings. 

This model for dark matter was originally introduced in refs. [15, 16, 22, 28]. Here we 
modify the underlying implementation of the model and significantly extend the analysis of 
ref. [16] by inclusion of all latest observational constraints. We also study more closely some 
benchmark scenarios for which we perform a new MCMC scan of the parameter space. We 
also analyze new consistency constraints arising from LHC-bounds on invisible and radiative 
Higgs decays. We find that the model can provide a naturally stable DM particle with a 
mass in the range moM ~ 30 — 800 GeV. Most of the allowed parameter space is within the 
reach of the next generation of direct and indirect DM search experiments. Also electroweak 
precision data, especially the Peskin-Takeuchi ^-parameter, provide stringent bounds on the 
model. 

We also show that the model gives rise to perfect 1-loop unification of all gauge couplings, 
including also the new technicolor interaction. The common unihcation scale and coupling 
are accurately determined: Mu ~ 2.2 x 10^^ GeV and a\j ~ 0.0304. The unihcation of 
the technicolor coupling is very sensitive to, and hence accurately determines the TC 1-loop 
IR-pole: « 340 GeV. Based on the QCD-analogue, this value strongly supports the 

natural TC scale Atc ~ 3 TeV. This result provides an important check on the internal 
consistency of the model. 

The paper is organized as follows: We briehy review the model in section 2 complete with 
determination of all gauge- and scalar couplings and the mass terms and a full implementation 
in the TC context. In section 3 we study the unihcation of the standard model and the 
technicolor gauge couplings. The dark matter analysis together with imposing laboratory 
constraints is described in section 4. Here we also show the results MCMC-scan of the 
model parameter space and show that our WIMP is compatible with all current observational 
constraints. In section 5 we consider the LHC constraints on higgs decays and radiative 
corrections to its mass. Finally in 6 we give our conclusions and outlook. 
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2 General model of dark matter 


We consider a model where the dark matter candidate is a mixture between electroweak 
singlet and neutral components of electroweak doublet and triplet fermions. Such candidate 
arises from a low energy Lagrangian of the form 

= Af,g + -^Ad.g + Af,H + ^Ad,H + , (2.1) 

where vC 4 f,g describes the gauge sector of heavy 4^^ lepton family and £Ad,g that of the SU(2) 
adjoint and singlet Weyl fermions and and /lAd,H detail their interactions with scalar 

fields. Finally £sm describes Standard Model, where Higgs may be either a fundamental or 
an effective doublet field. We will now briefly introduce the relevant terms in the Lagrangian 
following ref. [16]. 

2.1 Gauge interactions 

We denote the left handed heavy doublet by Ll = and the charged right handed 

singlet by Eji. Using the SM-like hypercharge assignments the kinetic and gauge interaction 
terms for these fields become equal to the corresponding terms in the SM: 

£4f,g = iLl ^ Ll + iEn ^ Er + £w + + l^A , (2.2) 

where the gauge currents are given by 

£w = ^ {W-Eli^Nl + W^Nli^El) , 

Cz = + {2sl - 1 )El7>^El + 2sIErj>^ErJ , 

Ca = -eA^E-f>^E , (2.3) 

where g is the weak coupling constant and Cw = cos6w and Sw = sin^vV) where 6w is the 
Weinberg angle. 

The Lagrangian for the left handed SU(2) adjoint triplet uj = and the right 

handed singlet Weyl fermions is 

^Ad,g = ioJ^a^D^u) + . (2.4) 

Here = (1 ,(t) and where cj* are the usual Pauli matrices, and the covariant 

derivative is We can go to the 4-component notation by defining 

charged Dirac spinors w~j^ = {w~ (rc^)^")^ and w~^ = {w'^ (r(;“)^")^ and neutral Majorana 
spinors wm = (tc^)!")"^ and /3 m = {Pa /3^“)^, where = {w^ =F iw‘^)j^/2 are the two- 

component charge eigenstates. In the 4-component notation the Lagrangian (2.4) becomes: 

^Ad,g = iWD pWD + ^WM pWM + ^Pm $ Pm 

+ g (VF+ wmYwd + w~WDl^WM - W^wdI^wd) , (2.5) 

where wr = w~jj. Note that the neutral field wm does not couple to neutral gauge boson, 
and that the singlet field Pm has no gauge couplings. This interaction pattern has important 
consequences for the observability of the DM particle in the model. 
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2.2 Scalar interactions and masses 

The masses and mixings of the new fields introduced above are determined by their couplings 
with the scalar sector, which we represent by an effective SM-like Higgs doublet H. First, 
the gauge-invariant interactions between H and the 4*^ family leptons and the neutral singlet 
are, up to dimension five operators, given by: 

Af,H = VeLlHEr + ypLiMf^R + ^{LlH){H^LI) + h.c., (2.6) 

where H = and ue, yp and Ann are dimensionless coupling constants and A ':$> v 

is yet some unknown scale related to the UV complete theory generating the full flavor 
structure of the model. After symmetry breaking the first Yukawa term generates a Dirac 
mass ruE = yE^lV^ for the charged lepton E, where v is the vacuum expectation value of the 
neutral composite Higgs held h. Second Yukawa term and the non-renormalizable dimension 
hve operator produce Dirac and Majorana mass terms for the neutrino mass matrix. 

Note that all interactions in Eqs. (2.2), (2.5) and (2.6) are invariant under Z 2 symmetry 
transformation, where 

E ^—E, N ^—N, /3 ^ —/3 and w ^—w. (2.7) 

We take this to be an exact symmetry of the model, which then forbids any couplings between 
the SM particles and the doublet or triplet helds. This symmetry is crucial for the stability of 
the dark matter in the model. We next allow all couplings between H and the SU(2) adjoint 
helds, which are consistent with the Z 2 -symmetry, again up to dimension hve operators: 

CAd,n = VwH^ujLl + + h.c., (2.8) 

where ui = and t“ = cr“/2 in terms of the Pauli matrices and the scale A is the 

same we introduced in Eq. (2.6). Finally, we include a gauge- and Z 2 symmetric interaction 
Lagrangian to provide a mass to the singlet held /3 

= ynSPP + h.c. (2.9) 

After symmetry breaking the interactions (2.6) and (2.8-2.9) give rise to a 3 x 3 mass matrix 
for the neutral Majorana particles N, wm and Pm (we drop the index M from wm and Pm 
here): 

-y / A^nn nT-Ntu nT-N/3 \ / AIr \ 

^mass = 2 {A!r, wr, Pr) TTT-n^ + h.c. (2.10) 

\ m^p m^p Mpp ) \Pl j 

where A/nn = Ann^^/A, = A,„„u^/4A, Mpp = V^yRVs, ruN^ = ywv/2y/2, m^p = y^vj^pi 
and vn^p = X^pv"^pA, where Vg is the VEV of the singlet held S. The lightest mass eigenstate 
of this mass matrix is stable by the Z 2 symmetry, and is thus identihed as the DM particle. 
Finally note that the Majorana mass is simultaneously the mass of the charged adjoint 
held Wr. 

^We assume that the bare mass term for /3 is zero even in the absence of a protecting symmetry principle. 
Whether such symmetry exists is related to a broader issue concerning the underlying flavor physics yielding 
also the effective couplings between the scalars and fermions. This is an interesting topic which, however, is 
beyond the scope of the analysis carried out in this paper. 


_ 4 „ 





2.3 Rotation to the mass eigenbasis 

The symmetric mass matrix appearing in Eq. (2.10) can be diagonalized by a unitary trans¬ 
formation U'^MU = m, such that the mass eigenvalues are rrii > 0. Using the notation 
Ql = {Nl, w^, /3l)^ Eq. (2.10) can then be written in the form 

^mass = , (2.11) 


where m is the diagonal mass matrix with positive mass eigenvalues. The corresponding 
mass eigenstates are Majorana fields given hy x = XL + XR = + U'^Q.r . This relation 

can be immediately inverted to give Or = Uxl and Or = U*xr- Using these relations 
and Eqs. (2.3) and (2.5) we then hnd the weak currents of the heavy leptons and the SU(2) 
adjoint fermions in the mass eigenbasis: 


Cf, = -j^W- Y, UuELl^^XrL + h.C. , 

^ i 



^Ad = Y ^D{Vi + lAi'y^'j'yf^Xi + h.c., 
i 


where 


( 2 . 12 ) 

(2.13) 

(2.14) 


U,- = 9(C/r,C/i,), Aij = Vi = ^{U 2 i) and A, = Q(U 2 i) (2.15) 

and Uij are the elements of the diagonalizing matrix U. Similarly, from Eqs. (2.6), (2.8) and 
(2.9) we can find the Higgs interactions in the mass eigenbasis: 

£-.x ^-^T,XiiSii + Pifhx, - 0 -E+ aiA)x. + ... . (2.16) 

i<j ^ i 

Here dots refer to terms which do not affect tree level matrix element calculations and mw 
is -boson mass. The various mixing angle and mass dependent coefficients are defined as 


Sij — 171]^^Aij “b (^Sij ‘2^M^^Dij TTliSij , 

Pij — Lfi-^piVij i{6ij ‘2'jMppEij , 

Sii — TTl-f^pAii 2MppDii -\- -^Tfli , 

Pii = —m^piVa + ^MppEii , (2-17) 

where m* is the i’th mass eigenvalue. The projection factors Uj and Aij are as defined in 
Eq. (2.15) and 

Dij = ^{UsiUsj) and Eij = ’AiUsiU^j). (2.18) 

Equations (2.12-2.18) contain all information needed to calculate the WIMP interaction rates 
relevant for the relic density and direct detection analyses. 
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2.4 Renormalizable and anomaly free implementation 

The dark matter model discussed above is not a consistent extension of the SM on its own 
as it suffers from quantum anomalies. However, the low energy Lagrangian (2.1) can be em¬ 
bedded e.g. into the context of a renormalizable, anomaly free TC model where electroweak 
symmetry is broken dynamically. One possible realisation is to take this to be the minimal 
walking TC, and this possibility was explored in ref. [16] (see also [15, 22, 28]). Here, to 
illustrate different possibilities, we consider an alternative realization. The complete list of 
new fields and their quantum number assignments is shown in table 1. In addition to the 



SU(3), 

SU( 2 )i 

U(l)y 

SU(iVTc) 

^2 

Ll 

1 

2 

- 1/2 

1 

-1 

E'k 

1 

1 

1 

1 

-1 

iO 

1 

adj. 

0 

1 

-1 

13 

1 

1 

0 

1 

-1 

9 

adj. 

1 

0 

1 

-1 

Ql 

1 

2 

1/6 

3 

1 


1 

1 

-2/3 

3 

1 


1 

1 

1/3 

3 

1 

m 

1 

1 

0 

3 

-1 

m 

1 

1 

0 

3 

-1 

G 

1 

1 

0 

adj. 

-1 


Table 1. The table shows the new states added to SM, and their charge assignments under the SM 
gauge group and the technicolor gauge group which we will consider to be SU(3). Also shown is the 
discrete matter parity which is even for the SM matter fields. 

fields introduced earlier (L^, oj and /3), the table shows the new strongly coupled sector 
responsible for the electroweak symmetry breaking. Of the new fermion fields, the techni- 
quarks Ql-, H)), rji and 772 , are gauged under a new vectorial gauge interaction SU(A^tc)- 
We set A^tc= 3. These elementary fermions form composite fields similar to the mesons and 
hadrons in QCD. We assume that only one doublet of the technifermions is gauged under the 
electroweak symmetry, while the remaining two Weyl fermions (771 and 772) are singlet under 
all SM charges. We assume that similarly to other SM-singlet fermions also 771 and 772 are odd 
under the ’’matter parity” Z 2 . Finally, there are two Weyl fermions, g and G, transforming 
in the adjoint representation of SU(3) of QCD and TC, respectively^. Both of these fields 
are assumed to be heavy and decoupled from low energy particle spectrum. These field are 
relevant neither for the dynamical symmetry breaking nor for the dark matter. However, as 
we will discuss in the next section, when g and G are included, along with the SU(2)-adjoint 
field w, the model also gives rise to excellent gauge coupling unification. 

At high energies the technicolor sector is described by the Lagrangian 

c^c = - QLiJ/>LQL - UrUPjPr - DrUP^Dr -gilpg, (2.19) 

where F^y is the field strength of the technicolor gauge field and Q = {U,D)'^. The covariant 
derivative D contains only the technicolor gauge field while the covariant derivatives Dl^r 
contain also the electroweak gauge fields. At low energies the strong dynamics is described 

^For concreteness we have assigned the value Z 2 = — 1 for these fields. 
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by an effective Lagrangian for composite mesons. Due to the different parities of the 
techniquarks, the low energy composites are S ~ QQ and a ~ r/r/. The former is the effective 
Higgs doublet which in our case is a composite field, and hence the model does not suffer 
from hierarchy problem. The field cr ~ S' + mg is another composite complex scalar, singlet 
under all SM charges. 

The low energy effective Lagrangian is 

£TC,eff = TVD^StD^S + - H(S, a), (2.20) 

where S = ((^ + iff • (t)/ 2 is charged under the electroweak interactions (ui are the Pauli 
matrices) and 

P(S, a) = mhrM^M + (2.21) 

where M = S 0 a. The real part, S, is identihed with the field introduced in Eq. (2.9) 
to provide mass to the fermion field /3, while for the pseudoscalar component tTs we have 
included an explicit mass term in Eq. (2.21). Even if the field tt* is light with respect to the 
intrinsic scale Axc; h can be heavy with respect to the masses in the dark matter sector. If 
fis was comparable to, or smaller than the fermionic DM mass, the dark sector of the model 
would be more complicated, containing both fermionic and bosonic components. While this 
is an interesting possiblity, we shall here assume that is heavy and consequently the DM 
is purely fermionic. A mass term for vr^ is expected to arise from the flavor physics providing 
the masses of SM matter fermions, and their origin can be accounted for by gauge dynamics 
(as in extended TC models). 

As a final remark, we note that with the above particle content the TC sector might be 
close to conformality [30]. However, coupling between the TC sector and the SM fields via 
the gauge and Yukawa interactions will move the theory away from the conformal window 
[31-33]. Therefore, we will assume the properties of the TC sector to be QCD-like. 

3 Unification 

3.1 Unification of the SM coupling constants 

Eor completeness we will first briefly review the argument for the unification of the SM 
coupling constants. More details can be found in ref. [15]. At one-loop the coupling constant 
an of an SU(n) gauge theory is given by 

a“^(^) = a-^{Mz) - ■ (3-1) 

The beta function coefficient is: 

bn = lT{R)N^f + ^T{R')N,b - yC 2 (G), (3.2) 

where T{R) and and T{R') are the Casimirs of the representation R for Nyjf Weyl fermions 
and of the representation R' for Neb complex scalars and C 2 {G) is the quadratic Casimir 
of the adjoint representation of the gauge group. Eor SM we have three coupling constants 
corresponding to n = 3,2,1. Requiring that SM coupling constants unify means that the 
three couplings are all equal at some scale M\j: a 3 (Mu) = a 2 (Mu) = ai(Mu) with ai = 
ajii? cos^ 6w) and a 2 = aj sin^ 6w, where c is a normalization constant that depends on the 
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choice of the unifying group. Here we shall use c = y^S/S, corresponding to the SM matter 
unified into SU(5). 

Using Eq. (3.1) we can now derive the following relation: 


B = 


bs - h 
b2 - bi 


aja-i — sin^ 0w 

(1 + c^) sin^ 0w — 


= 0.721 ±0.004, 


(3.3) 


where the Weinberg angle 9w and weak and strong coupling constants were evaluated at 
the Z-mass scale, using values from ref. [34]: sirP' 9w{Mz) = 0.23126 ± 0.00005, a~^{Mz) = 
127.940 ± 0.014, a 3 (Mz) = 0.1193 ± 0.0016 and Mz = 91.1876 ± 0.0021 GeV. 

The hypercharge assignment of our model renders the Technicolor sector identical to 
one extra SM generation from the electroweak interaction viewpoint. In addition we have one 
strongly interacting adjoint Weyl fermion, which affects the running of the QCD coupling 
and one weak triplet affecting the running of a 2 - The group factors are T{R) = 1/2 for 
the fundamental representation, T{G) = n for the adjoint SU(N)-representation and T{R) = 
g2y2 _ (3/5jy2 U{ 1 )y hypercharge gauge group. With Ng generations of ordinary 

fermions we then find: 


bi 


h2 


bs 


4.. 4 

-Na ± - 
3 ^ 3 

4.. 14 

-Na - 

3 ^ 3 

4 


6f“ + - 

^ 30 

bf' + ^ 


lSM 


± 2 . 


(3.4) 


Note that the differences bi — bj are independent of Ng, because they can not be affected by 
states forming complete representations of the unifying gauge group [35] . It is now clear that 
the SM does not unify since ~ 0.53. However in our model ~ 0.722, which is 

generously within one sigma of the extremely tight constraint (3.3). In fact the unification 
mass and coupling are very precisely determined by the coupling constant unification. If we 
define a chi-squared function 


X^(Mu,au) 


E 


2=1 


Aaf 


(3.5) 


where Aa* are the observational errors for each coupling given above in the text (these errors 
propagate essentially as such to the unification scale), we find that at la-level: 

Mu = (2.20 ± 0.03) X 10^® GeV, and au = 0.03042 ± 0.00002 . (3.6) 


In any grand unified theory nucleons are expected to decay via the exchange of gauge 
bosons with GUT scale masses. Schematically the partial decay width of the proton into a 
generic channel containing a meson and a lepton is 


r = 7QCD7GUT, (3.7) 

where 7gut contains the details of the underlying unified theory and 7qcd contains the QCD 
parameters and the effective low energy constants parametrizing the hadronic matrix element 
relevant for the decay in question. For example, for a simple decay mode via a massive gauge 
boson exchange, assuming M « Mu, 7gut ~ Ofy/M^. The precise details of course depend 
on the particular GUT model. We do not pursue such model building here; see e.g. [36, 37]. 
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Figure 1. Left; the running of the gauge couplings in the SM. Note the logarithmic scale. Left: the 
running of all four couplings in the MWTC-DM model under consideration, including the TC-coupling 
Q! 4 . The ellipses show the 1- and 2-a contours on unification scale and the unified coupling derived 
from the x^-distribution (3.5). For the TC-coupling 04 we took 331 GeV < < 351 GeV. Note 

the linear scale on the latter plot. 


To obtain a parametric estimate, consider 


tn 


1 n 

"n 9 95 

i niN 


(3.8) 


where f-^ = 0.131 GeV, ttin is the mass of the nucleon and the hadronic low energy 
constant. This must be determined from the lattice [38] and is subject to relatively large 
uncertainties [39, 40]; for p —)• the estimates for the value of oat range from 0.003 GeV^ 

to 0.03 GeV^ [39]. Using the value 0.01 GeV^ compatible with the lattice calculation [39] 
and numbers from Eq. (3.6) results in Tjy ~ 10^^ y, which is compatible with the current 
bound from the Super-Kamiokande tn > 10^“^ y [34]. 

We plot the running couplings in figure 1 for the SM (left panel) and for the current 
model (right). The latter plot was created with linear M\j scale and zoomed to the unifi¬ 
cation coupling to reveal the almost perfect one-loop unification in our model. The tight 
error bars on the unification mass and the unified coupling can be used to make a formally 
accurate prediction for the value of the QCD-coupling at electroweak scale. Running the 
QGD-coupling backwards from the unification scale gives: 


asiMz) = 0.1120 ±0.0003, 


which is consistent with but much tighter than the current observational limits^. 


(3.9) 


3.2 Unification of all couplings including 04 = otc 

Above we considered only the unification of the SM coupling constants. In our model we 
have an additional gauge coupling related to the strong Technicolor interactions and it would 

^Note that the running of the SM couplings will be affected by the strongly coupled TC sector at the scales 
below C>(TeV): the composite spectrum of technihadrons charged under the electroweak interactions will feed 
into the evolution and may affect the precision of the above result. 
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be more satisfying to have a unification of all four coupling constants. We now show that 
this indeed quite naturally takes place in our model. ^ 

With the current particle content, shown in Table 1 we have: 


64 = (4 + 2) X ^ + 2-11 = -7. 

O 

Let us now require that 0:2 and 0:4 unify at M\j. This implies that: 

sin^ 9w — 0/0:4 ^4 — &2 


C = 


— (1 + c^) sin^ 9w b2 — bi 


Theoretically, the particle content of the model gives: 

64 — 62 


C'th = 


62 - 


1.0556. 


(3.10) 


(3.11) 


(3.12) 


The “experimental” value of C depends on defined as the scale at which the 

inverse 1-loop coupling vanishes: 04 ^(A;^’/,*°°^) = 0. Using = Mz gives « 

1.0055 ± 0.0006, which is formally about 90 sigma’s away from the theoretical value. Since 
< Uth and 64 is negative, the problem is that 04 undershoots the unification value. 
We can improve the situation either by adding new fermion or boson helds to the particle 
spectrum to achieve a slower running, or by use of a larger A;^’/,*°°^ to get an effectively 
positive a^^{Mz) initially in Eq. (3.11). Let us consider the latter option. From Eqs.(3.11- 
3.12) it is easy to see that the scale needed to achieve unification between 0:4 and 0:2 

is given by: 

, 1—loop 
^TC 


064 , An 




Mz 


= [Cth - {<? - (1 + <?) sin^ Qw) ■ 


(3.13) 


Using the experimental input values for couplings and for the above computed values for Cth 
and propagating all errors throughout, we find: 


A, 


1—loop _ 


TC 


= 341 + 5 GeV. 


(3.14) 


The right panel of figure 1 illustrates how the complete 4-couphng constant unification takes 
place at the scale and the coupling given by Eq. (3.6), when A^’/!°°'’ is chosen according to 
Eq. (3.14). Of course the energy scale where TC dynamics becomes strong is above the scale 
where the 1-loop coupling diverges. For example, for QCD itself, the simple 1-loop 
running Eq. (3.13), yields « 57 MeV, while the typical scale for QCD is about an 

order of magnitude higher: Aqcd ~ 700 MeV. If this holds also for TC, then the unification 
condition Eq. (3.14) predicts that 

Atc ~ 3 TeV . (3.15) 

This agrees very well with what one would naively expect: if we write Atc ~ dvrEn-; where 
En- is the technipion decay constant, then the unification condition (3.15) sets E^r ~ 250 GeV. 

Even though such scalings are simply naive dimensional analysis, it is encouraging 
that unification of all couplings is not only possible in the model, but also determines the 

"^The results are dependent on the normalization of the hypercharge, i.e. the factor c. In principle its value 
is determined by the particle content and the structure of the unifying algebra; we assume the value c = \/3/5 
throughout here. 
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dynamical symmetry breaking to occur at the TeV scale. Even if our analysis is complete at 
1 -loop level, the results must be taken with a grain of salt: detailed results may be modified 
by higher order corrections and by threshold effects. Moreover, we have not addressed the 
underling nature of flavor dynamic which gives rise to the couplings between the effective 
scalars and fermion fields, as this dynamics is decoupled from the other gauge interactions 
at 1-loop level. As an initial exploration of this model, our results nevertheless provide an 
interesting benchmark scenario for more detailed investigations in the future. 

4 Experimental constraints on low energy theory 

At low energies our model is parametrized by seven dimensionless couplings ys, yp, yji, 
Ann) ^-wp sud thrcc scales v, Vs and A. From these one can easily work out the entries in 
the mass matrix Eq. (2.10) of neutral fields and the mass of the new charged state E. From 
naive dimensional analysis we infer that y* < dvr and Xij < (dvr)^. Furthermore, assuming 
Vs ^ V = 2d6 GeV and A ~ O(TeV), we find it reasonable to adopt the following prior ranges 
for the Lagrangian masses: 

\Nkj\ < 3000 GeV; \mij\ < 2000 GeV and 200GeV < niE < 2000 GeV . (d.l) 

We scanned this parameter range using Monte Garlo Markov Ghain (MGMG) methods [16]. 
For each set of input parameters we diagonalize the DM mass matrix numerically, find the 
mass eigenvalues rrii and the diagonalizing matrix Uij and identify the lightest eigenstate as 
the WIMP. We then check that the WIMP is stable, i.e. that it is the lightest of all states 
transforming nontrivially under the Z 2 symmetry. 

The data shown in our result figures is compatible with the experimental and obser¬ 
vational constraints from oblique electroweak precision data, the Z-boson and Higgs bo¬ 
son invisible decay width limits, cross section constraints from DM direct detection LUX, 
XENONIOO and PIGO experiments as well as DM indirect detection constraints from Ice- 
Gube, Super-Kamiokande and FERMI-LAT telescopes and from the AMS-02 experiment. 
For the data passing these tests, the DM relic density is calculated numerically and checked 
to be consistent with the most recent observations [41]. 

We do not require that our model provides the total observed abundance of the DM, 
inferred from the most recent GMB observations: = 0.1193 (±0.0014) [41]. Instead, 

we impose this as an upper bound and compute how large a fraction of the total DM-density 
each parameter set is able to produce, defined as frei = ■ We accept models also 

with subleading DM in the interval: 


0.05 </rel < 1.01. (4.2) 

This criterion affects the direct and indirect DM search constraints on WIMP-nucleon cross 
sections, which usually are given assuming that /rei = 1- However, as long as different 
DM-components are weakly interacting, they all cluster roughly the same way, and a given 
subleading DM should make up only a fraction /^ei of the DM density in all cosmological 
substructures. We then constrain such subleading WIMPs using a scaled effective cross 
section [8, 9]: 

^SD,SI = /rel CrsD.SI < O-fand , (4:.3) 

where usd refers to spin-dependent and usi to spin-independent channel and Ubnd is the 
bound from a given experiment. We imposed direct search bounds from XENONIOO [42], 
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LUX [43] and PICO [44], as well as the indirect search bounds from IceCube [45, 46] and 
Super-Kamiokande [47, 48]. For explicit expressions for cross sections and for more detailed 
discussion of the implementation of these constraints see [16]®. Both SI- and SD-constraints 
are relevant for our model in different regions of the parameter space. However, the SI- 
constraint is typically the stronger one. There are particular cases where our DM particle 
has essentially but a dominantly pseudo-scalar coupling to the Higgs boson. In such case 
the WIMP-nucleus interaction is momentum transfer dependent and strongly suppressed. 
These solutions may avoid detection by any of the DM search programs currently under 
construction. 

The new doublet and adjoint SU(2) states in our model, as well as the new states in the 
TC sector are charged under SU(2) and hence contribute to oblique S and T-parameters [49]. 
Explicit expressions of these contributions can be found in the appendix of ref. [16]. Here we 
use the experimental constraints [34]: 

S = 0.00 ± 0.08, and r = 0.05 ±0.07. (4.4) 


which include a 90% correlation between S and T as given by [34]. There are many other 
bounds coming from collider experiments. First, there is a direct LEPH-bound on any 
charged particle i coupling to X-boson: m* > 104.5 GeV. LHC mass limits, while not as 
straightforward to implement, are typically much stronger. In our analysis we have used 
conservative bounds 

iTiojo > 500GeV. (4.5) 

The X-boson invisible decay width imposes a constraint on any particle with m < Mzl‘2‘- 
The current bound from LEPH is r(X —)■ inv.) = (2.984 ± 0.008)r(X uu) [50]. As the best 
fit value is already 2cj below the SM prediction, we allow at most one standard deviation 
from new physics, which implies a bound 

/ Am ^ \ 

6z^\Uu\Ul -< 0.008. (4.6) 

' 7Tt rj ' 


This rules out any WIMP with tudm < mzl‘^ and a significant Ni component. Furthermore, 
if the WIMP is lighter than then also Higgs could decay to a pair of WIMPs. The 

invisible Higgs branching fraction Ri is constrained to be [51-53]: 


Ri = 


Th, 


,DM 


Fh,DM ± PsM.tot 


<0.17, 


(4.7) 


where FsM,tot is the total Higgs decay width in the SM and Fh,dm = (Girm^/SV27r)( 
l-PiiPA); where /3j = (1 — 4mf/m^)^/^ and the index i refers to the WIMP as the lightest 
of the mass eigenstates. The bound (4.7) assumes SM-like Higgs-gauge and Higgs-fermion 
couplings. It would be relaxed to Ri < 0.26, if one allows Higgs and SM gauge fields to have 
non-SM-like couplings to photons and gluons [53]. However, in our analysis we always use 
the stronger constraint (4.7). 

®Here we have improved our analysis related to IceCube and Super-Kamiokande limits by taking the DM 
annihilation branching fractions to different channels into account when imposing the constraints. We use 
W'^W~ limits for annihilation channels W^W~,ZZ and Z/i, and and bb limits as they are, including 

proper branching fractions in all channels. 
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Figure 2. Scatter plot of the models passing all existing constraints as a function of the DM-mass 
and the predicted WIMP-nucleon Sl-cross section. Colors represent the value of /rei for each model 
as indicated by the bar on the right. Also shown are the current XENONIOO [42], and LUX [43] 
limits as well as the predicted reach of the XENONIT experiment [54]. Large red dots show our five 
benchmark models. 


Indirect observations are also sensitive on WIMPs annihilating in Galaxy center and in 
Galaxy halo. However, current limits on signals from neutrino detectors [55-57] are not strin¬ 
gent enough to be of use here. These annihilations could also create potentially observable flux 
of gamma-rays in the FERMI-LAT data [58-62]. However, our light WIMPs (mnivi < fnw) 
tend to annihilate to lighter fermions in a velocity suppressed p-wave via Z-boson resonance 
giving only a very weak gamma-ray signal. Thus, in our model the annihilation cross sec¬ 
tions are below the current FERMI-LAT limits [58-62] in the most constrained low DM mass 
region. Only the FERMI-LAT limits from the Milky Way dwarf spheroidal satellite galaxies 
(dSphs) [63] are sensitive enough to constrain a small part of our model parameter space. 
In WIMP mass region mw ^ w-dm ^ 100 GeV these limits [63] cut away a few points. For 
^DM ^ 100 GeV these limits have no impact. In all our results we show only data which 
passes these FERMI-LAT constraints. Finally, even though the annihilation cross section for 
subdominant DM in dSphs is larger than the canonical thermal relic cross section, the DM 
density is also expected to be smaller, scaled down by factor /rei- This suppresses the flux of 
gamma-rays by a factor making these constraints for subdominant DM even milder than 
for canonical thermal DM. 

Finally, we imposed the constraints for DM annihilation cross section in the bb channel 
from ref. [64] (left panel Fig. 4), derived from the new AMS-02 [65] data on the secondary 
astrophysical antiproton to proton ratio. These limits cut away the few otherwise remaining 
points in the mass range todm ~ 65 — 80 GeV. 
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Figure 3. Scatter plot of the models passing all existing constraints as a function of the DM-mass and 
the WIMP-proton SD-cross section. Shown are also the best current constraints direct [44] (PICO) 
and indirect searches [45, 46] (ICECUBE) and [47, 48](Super-Kamiokande). In models indicated with 
dots WIMPs would have reached equilibrium between capture and annihilations in the sun and the 
models indicated with crosses they have not. For more detailed treatment of the indirect observation 
channels see [16]. 


4.1 Results of a generic MCMC scan 

In figure 2 we show the distribution of the models that passed all tests in our MCMC runs 
as a function of the DM-mass, the effective WIMP-nucleon Sl-cross section and the relative 
relic abundance /rei, whose value is indicated by the vertical bar to the right of the plot. The 
advantage of our using the effective cross section here is that one immediately sees how much 
a given direct search experiment needs to improve its sensitivity in order to rule out a given 
set of parameters. It is still easy to find acceptable models, in particular with a subleading 
DM. However, most of the allowed parameter space, including all our benchmark models to 
be defined below, is within the reach of the next round of the direct search experiments, 
which improve the current bound on Ug® by a factor of ~ 50. 

Let us now very briefly explain the data in figure 2. When tudm ^ 80 GeV, the DM 
particles annihilate into light SM fermions. A vertical cluster of points around tudm ~ 45 
GeV corresponds to the Z-boson resonance. Another, weaker cluster, around tudm ~ 60 GeV 
(containing only blue points) corresponds to the Higgs boson resonance. Solutions at range 
moM ~ 65 — 80 GeV are exlcuded by AMS-02 constraints as was previously explained. For 
^ 80 GeV, new annihilation channels, first to W~W~^ and then to ZZ, Zh, hh and tt 
open up sequentially and begin to dominate the cross section. None of these channels have 
equally striking signature as do the resonances. 

In figure 3 we show the projection of accepted points as a function of the spin de¬ 
pendent interaction cross section cjgp p and the WIMP mass. It is evident that the exist¬ 
ing Sl-constraints are stronger than the SD-constraints. However, the IceGube and Super- 
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Figure 4. Scatter plot of the models passing all existing constraints in the plane of precision elec- 
troweak parameters S and T. The ellipses are the experimental Icr, 1.6a and 2.6a confidence contours 
for S and T [34]. Note that the benchmark model 2 is as good a fit to the data as is the SM. 



Mnn 


Mi3i3 



-^ui/3 

ruE 

Ml 

3.12-il.81 

15.7-i0.37 

0.60-bi0.04 

0.13-bi0.72 

-0.89-i0.ll 

-0.45-bi0.07 

7 

M2 

2.18-i0.17 

8.92-bi0.14 

0.92-bi0.04 

-0.02-i0.38 

-0.29-i0.13 

-0.69-bi0.23 

6.5 

M3 

7.33-i0.85 

6.59-bi0.49 

1.42-bi0.03 

0.51-bi3.05 

0.52-i0.41 

-0.70-i0.58 

12 

M4 

9.39-i0.74 

7.92-i0.06 

2.70-bi0.14 

-0.21-bi2.84 

-0.88-i0.73 

-0.89-i0.53 

18 

M5 

13.1-bi0.14 

12.7-bi0.25 

8.12-i0.06 

-0.16-i0.13 

-0.38-i0.ll 

0.38-bi0.03 

15 


Table 2. Shown are the inpnt mass parameters of the benchmark models M1-M5 along with the 
mass of the new charged doublet state ruE. All masses are given in nnits 100 GeV. Note that the 
mass of the new adjoint state ojd equals with the Majorana mass 


Kamiokande limits, that concern the WIMPs that accnmnlate in the core of the Snn and 
then annihilate to kP-bosons, have some constraining power in the mass region ttidm ^ 200 
GeV. At first sight it seems that the latest Super-Kamiokande results, concerning WIMPs 
that annihilate into tan leptons, also have some constraining power. Note however, that the 
points above the Snper-Kamiokande T'’'r“-hne in Fig. 2 are not excluded. This is because 
the Super-Kamiokande constraint assumes that WIMPs annihilate into tans with a branch¬ 
ing ratio Br,-,- = 1. Here in general, and in the accepted models falling above the Super- 
Kamiokande constraint in particular, the branching ratio to tans is much less than one. For 
DM masses below W^W~ threshold our WIMPs annihilate dominantly to 66-channel, and 
the T^T~ branching is only ~ 0.05 (for branching ratios in specific benchmark models see 
table 3). Furthermore, this ratio only decreases once other annihilation channels open for 
iTT-DM A mw- Due to these small annihilation branching fractions, the Super-Kamiokande 
limits are currently very little constraining. 
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Figure 5. Left: shown is the result of a restricted MCMC-scan in the complex mass parameter muip 
around the benchmark point Ml in Table 2. Light gray dots show the projections of the points to the 
complex mass plane. Right: the projection of these points in (todm, cTgf )-plane. 


4.2 Benchmark models 

A generic MCMC scan over the entire prior range gives a good idea of the constraining power 
of the different observations. However, because of the high-dimensionality of the parameter 
space, such scans do not reveal the hirer details of how acceptable models are distributed. 
In particular it appears that there are but a few sets of parameters that give /rei ~ 1- For 
this reason we selected hve benchmark points from the accepted MCMC data sets and made 
new runs with restricted priors in their neigborhoods. The selected models are labelled as 
M1-M5 and shown by large red dots in the scatter plots 2-6 and 8. The corresponding 
central parameter values are given in table 2. In table 3 we show the ensuing DM mass, 
relative DM-abundance, precision electroweak parameters S and T, the effective WIMP- 
nucleon cross section a^, the tau-branching ratio Br,-,-, the contribution to Z-width and 
the invisible Higgs decay fraction Rj for these models. Note that models 2-5 have large and 
positive S and T-parameters. This is a generic feature in our model, due to the fact that 
the precision variables, and S in particular, get a large positive contribution from technicolor 
fields (see table 1). As is evident from hgure 4, a strict bound on the 5-parameter 5 < 0.1 
could rule the model out completely. However, the contributions from flavor extensions are 
subtle and may quantitatively affect the analysis [66, 67]. We give no values for 5 and T 
for model 1, because our precision data analysis [16] is not applicable for ttidm < Mzl2. 
However, we expect that the strongest bound comes from the Z-decay width in this region. 

In the left panel of figure 5 we show a scan of parameters around the benchmark Model 
1. We hxed all parameters as given in table 2 except My^p, which was allowed to vary freely 
in an MCMC scan starting from the benchmark value. This scan reveals a continuous, but 
constrained domain of parameters giving a right, or very closely right DM abundance. In the 
right panel we show how these points are all tightly concentrated in the the DM mass-effective 
Sl-cross section diagram. This example shows that the apparent deficit of points with /rei ~ 1 
in Figs. 2-4 may give a too pessimistic idea of the density of acceptable models. Figure 6 
shows a result of an MCMC-run around the benchmark model 3, fixing all parameters but 
the complex mass Mnn- First, we observe that /rei is almost independent of Mnn- This is 
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moM 

/rel 

S 

T 

ett 



Ri 

Ml 

38 

1.00 

- 

- 

2.2 X 10-4*^ 

0.051 

0.0008 

0.12 

M2 

84 

0.96 

0.11 

0.13 

1.7 X 10"^® 

0.0012 

- 

- 

M3 

137 

1.01 

0.18 

0.14 

1.3 X 10"^® 

7.7 X 10"® 

- 

- 

M4 

246 

1.01 

0.16 

0.15 

9.1 X 10-^® 

1.3 X 10-® 

- 

- 

M5 

806 

0.87 

0.22 

0.19 

3.7 X 10"^® 

9.9 X 10"^ 

- 

- 


Table 3. Shown are the values of DM mass todm, relative relic density /rei, precision S and T 
parameters, the effective cross section dg®, the contribution to Z-width 6z and the invisible Higgs 
decay fraction Ri for the benchmark models. A dash indicates that the bound is not relevant for the 
model in question. 




Figure 6. Left: shown is the result of a restricted MCMC-scan in the complex mass parameter ttinn 
around the benchmark point M3 in Table 2. Right: the projection of these points in (S', T)-plane. 


because DM is always mostly /3-like and so its mass and couplings do not depend much on 
the NN-entry of the mass matrix. Second, only a thin line of acceptable solutions are found. 
The reason for this is the T-parameter, which gets a large contribution from the doublet¬ 
like states, and the contribution from the A^-like neutral state must accurately cancel the 
contribution from the charged ill-state [22] . For a fixed m e this works only for a very narrow 
range in the mass of the A^-like state, which is essentially set by Mnn- This explains why 
points with /rei ~ 1 are relatively sparsely distributed in the generic MCMC plots: our full 
parameter space has many dimensions (thirteen) and good models are forced to lie on narrow 
low-dimensional strips, which are hard to locate in a full parameter space scan. However, 
when good solutions are found, one in general finds continuous sheets of acceptable solutions 
in their immediate neigborhood. 


5 Additional bounds: Higgs to 77 and a light composite scalar 

As already emphasized, our generic DM setting provides a very attractive solution for the 
particle dark matter problem, the hierarchy problem and the gauge unification. However, 
before concluding, we still need to discuss two additional constraints. First, the bound on 
higgs decay to two photons and second, the lightness of the higgs mass in the effective field 
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theory picture underlying the model. We find that both issues can be addressed within the 
strongly interacting sector of the model and do not have consequences for the DM physics. 
Even though these bounds are of no concern for our main results on dark matter, we discuss 
them for the benefit of the overall consistency of the model. A reader not interested in these 
details may skip directly to the conclusions. 

5.1 Higgs to 77 

There are a number of new particles in our model, which can influence the Higgs boson decay 
widths. As we have already discussed, if our DM is light enough, todm < it opens 

up a new invisible Higgs decay channel, which would change the predicted higgs branching 
ratios. In addition, all new charged particles contribute to Higgs decays to two photons via 
loop corrections. 

The contribution from new TC singlets, the new heavy electron and the new charged 
state nested in the SU(2) triplet, are easy to compute. The sector which is not singlet 
under TC is more model dependent. To estimate this effect we consider the contributions of 
heavy massive resonances within the sigma model-like effective theory describing the Higgs 
sector at low energies. The loop corrections induced by vectors to 77 -amphtude are much 
larger than those from scalars with a similar coupling strength, and so vector resonances may 
easily dominate the effective H 77 coupling. Concretely, we consider a simple setup with one 
extra charged vector resonance, described as a new effective massive IT'-boson. The relevant 
effective Lagrangian then is 

V ^ V 

_ y hW'-w'^'- . (5,1) 

y V ^ 

f 

The sum / runs over the charged SM fermions and the new charged techni-singlet 
fermions. In practice these are the top quark, the new heavy electron E and the states. 
The effective fermion and EW vector boson couplings of the Higgs are denoted by c/ and cw 
so that for SM we have Cf = cw = 1. In our model 1 and c^± = 2, while the couplings 
to the top quark and VE-bosons must be inferred from the LHC data. Generally, on the basis 
of extrapolating from QCD and explicit model calculations [67, 68 ], it is expected that in 
this type of a model cw = cz « 1 and cj « 1 which is also confirmed by fits to the LHC data 
[ 68 ].® 

For simplicity and for less model dependence we have written the Lagrangian in terms 
of mass eigenstates. We also neglected the mixing of VE-and lEhbosons, which in general 
could lead to a correlation between lE-and lE'-boson couplings. With these assumptions all 
non-conventional effects due to charged resonances are modelled by the coupling factor cy/' ■ 

The full 1-loop decay width of Higgs to two photons is 

= (5.2) 

®Note that the Higgs coupling to gluons is SM-like, since we assume that all the new particles coupling 
with the Higgs are color singlets. 
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where the effective coupling is given by 


QH'y'y 


a 

8tt 


'^CfNiQ‘jFi/2{Tf) + cwFi{tw) + cw'Fi{tw') 
f 


(5.3) 


The color factors of non-standard model fields are = 1 and the loop factors Fi^tj) 

for fermions (i = 1/2) and for vector bosons (i = 1) are standard, 

F'i/2{Tf) = - 2 t / [1 + (1 - t /)/( t /)] , 

Fiirw) = 2-1- 3tw + 3Tiy(2 - Tw)f{Tw), (5.4) 


where 


fi^j) 


arcsin' 


2 1 


log 




1 2 


— ITT 


if Tj > 1 
if Tj < 1 


(5.5) 


and Tj = 4mj/mjj. For tj larger than unity, the loop-factors quickly reach the asymptotic 
values Fi/ 2 (oo) = —4/3 and Ti(oo) = 7. Using the asymptotic value is a good approximation 
for all new fermions and even for the top quark the relative error of this approximation is 
about three per cent. Consequently, since the techni-resonances are expected to to be in 
O(TeV) mass range, the exact value of the W mass is not relevant and F\{m’\Yi) ~ 7. 

The effective coupling cw', can now be constrained by existing bounds on the ratio of 
the Higgs branching fractions to two photons i.e. BR.y.y/BR®^. We extracted a 2iT-limit 
for this quantity from the left panel of Fig. 5 of ref. [53]. The experimental bound is then 
converted to a constraint on cw' assuming no invisible decay channels. The result is shown 
in Fig. 7. The rather loose bound 0.52 < c^/ < 0.85 indicates that our model can be made 
consistent with the data by rather modest and reasonable assumptions about the structure 
of the effective low-energy theory. 

Another possibility would be to include the TC dynamics by considering the elementary 
fields, i.e. charged techniquarks, inside the loop and interacting with the TC-Higgs. We have 
checked, adapting the resuls of [68] to our SM-like hypercharge convention, that that this 
leads to quantitatively similar and consistent results with the effective model approach we 
have detailed here. 


5.2 Dynamical Higgs boson mass scale 

Then we discuss how a light composite Higgs relates to the intrinsic dynamical scale of the 
model. The traditional expectation of a heavy scalar resonance assumes TC dynamics in iso¬ 
lation. It is known that couplings with the EW gauge currents and extended flavour sectors, 
in particular the top quark, affect this conclusion [66, 67, 70]. In addition, our model contains 
several new heavy fields that couple to the Higgs and contribute to its vacuum polarization. 
We estimate this effect quantitatively by use of a perturbative one-loop approximation. As 
usual, we induce a cut-off A to evaluate the relevant 1-loop integrals. However, since A is 
not very large here, we keep nonzero masses for the heaviest particles in the loops. This 

^Recent results from Atlas collaboration [69] favor a slightly larger 77 -branching fraction, BRT..y/BR®l/ > 1. 
As is evident from Fig. 7, this situation is naturally realized in our model. 
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Cw' coupling factor 

Figure 7. Shown is the 2cr-limit for the coupling factor cw'- The allowed region is between the 
horizontal lines. The factor c^i is the effective coupling between techni-vector resonance and the 
techni-Higgs. 


procedure gives: 


+ [6m^ + 3m^ — 12mj] 


A2 






1 - 


A^ mf 


(5.6) 


Here mjp is the intrinsic dynamical TC Higgs boson mass, which we attempt to estimate. 
The cut-off dependent part in the first line includes corrections from the relevant SM particles 
and the sum in the second line the corrections from the new heavy fermions. These include 
the charged fermions E and co^ with factors Je = ^ and f^± = 4 and the two heaviest 
neutral Majorana fermions with a factor /^. = 4 for each®. We approximate the masses of 
the neutral states by M^n and M^^,, which is a reasonable approximation, as the Dirac mass 
terms are typically small compared to the diagonal Majorana masses in the original WIMP 
mass matrix. Thus the mixings are small and finally the contribution from the relatively 
light WIMP is suppressed compared to the other two Majorana states. 

The cut-off is provided by the TC-scale, which here, consistently with the unification: 
A ~ Atc ~ 3 TeV (see Sec. 3). For concreteness we used A = 1.5 TeV. We can now turn 
Eq. (5.6) around, setting tuh = 125 GeV and solve it for the dynamical Higgs mass rrifp 
for all models in our data set. The result is shown in Fig. 8. Models with 2mDM > 
give mjp ~ 1-2 TeV. This is of the right order of magnitude, since one would expect that 
in a conventional technicolor model mjp < dvru ~ few TeV. For models with tudm < Mzl2 
the small mixing approximation made in using Eq. (5.6) breaks down and the results for 
mjp cannot be trusted. Anyway, these arguments strongly suggest that the Higgs boson can 
indeed be naturally light in our model. 

®The factor /i = 4 follows from the expansion of the Higgs interaction term (1 + h/v)'^ —>■ 2hfv. This factor 
of two in the coupling gives the factor of four in the loop. 
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Figure 8. Scatter of the one-loop estimate for the intrinsic TC Higgs mass mjp for our accepted 
models. The colour coding again indicates the relative DM density. 


6 Conclusions 

We have analysed a model for dynamical symmetry breaking, dark matter and gauge coupling 
unification in light of the most recent observational and experimental data. The model 
setting is nontrivial, but very tightly constrained theoretically. Since the model does not 
contain any light fundamental scalars, there is no hierarchy problem. We have shown how 
the model gives rise to perfect 1-loop unification of all gauge couplings, including the new 
technicolor interaction, at a common unification scale M\j = 2.2 x 10^® GeV and the unified 
coupling a\] ~ 0.0304. Moreover, unification determines the scale of the 1-loop IR-pole of 
the TC-coupling ~ 340 GeV. By the QCD analogue, this is consistent with the naive 

expectation of the the TC interactions becoming strong around the scale Atc ~ O(TeV). 
Thus unification is not only possible, but it actually supports the existence of a TeV-scale 
strongly interacting dynamical sector. 

The essential part of the spectrum are the fermion fields transforming under adjoint 
representations of the gauge group. In terms of their quantum numbers these fermions are 
identical to the gauginos which arise in supersymmetric setting. Hence, it is natural to 
entertain the thought that the model is a low energy realization of a supersymmetric theory. 
The existence of a strongly coupled sector responsible for the electroweak symmetry breaking 
naturally decouples supersymmetry breaking from the electroweak physics, hence removing 
the little hierarchy problem [71-73]. Consequently the scalar superpartners can all be very 
heavy with masses around or above the unification scale Mu. More detailed model building 
and investigation of resulting phenomenology provide interesting further research prospects. 
It would be also desirable to have more detailed theory which explains the effective scalar- 
fermion couplings in our low energy lagrangian and the emergence of fermion mass patterns. 
The details of underlying flavor physics likely require an extensions of the technicolor gauge 
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dynamics. 

In this paper onr main analysis concerned the dark matter sector of the model. The 
dark matter candidate arises from mixing of three neutral fields: one gauge singlet, a neutral 
members of an SU(2)-doublet and an SU(2) triplet. The stability of the lightest of these 
fields is guaranteed by a discrete Z 2 -symmetry. The most essential parameters for the dark 
matter in model are the entries in the effective mixing mass matrix of these neutral states. 
We performed a generic MCMC scan of the model parameter space constraining the model 
by the most recent bounds following from the accelerators as well as direct and indirect 
DM-searches. 

We also introduced several explicit benchmark cases to illustrate typical features of vi¬ 
able models and performed limited range MCMC runs in the neighborhood of the benchmark 
points to study the allowed phase space in more detail. We found that there are large contin¬ 
uous regions of parameters for which the model can provide a naturally stable DM particle 
with a mass in the range tudm ~ 30 — 800 GeV. 

We conclude that the model is viable in light of existing data from collider experi¬ 
ments and cosmological and astrophysical observations. However, future experiments have 
excellent possibilities to probe the model further: Most of the available parameter space is 
within the reach of the next generation of DM search experiments. Also a significant shift 
of the observed precision electroweak parameters towards their SM-values, in particular of 
the Peskin-Takeuchi S'-parameter, could rule out the model as a source for DM. Finally, the 
experimental results on the proton decay will constrain the unification aspects of the model. 
Of course the interesting possibility is that one or the other of these observations would 
provide the first evidence of particles compatible with the low-energy spectrum predicted by 
the model. 
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